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Abstract 

Background: Slow-adapting type I (SA-I) afferents deliver sensory signals to the sonnatosensory cortex during 
low-frequency (or static) nnechanical stinnulation. It has been reported that the sonnatosensory projection fronn 
SA-I afferents is effective and reliable for object grasping and manipulation. Despite a large number of neuroimaging 
studies on cortical activation responding to tactile stimuli mediated by SA-I afferents, how sensory information of such 
tactile stimuli flows over the somatosensory cortex remains poorly understood. In this study, we investigated tactile 
information processing of pressure stimuli between the primary (SI) and secondary (Sll) somatosensory cortices by 
measuring effective connectivity using dynamic causal modeling (DCM). We applied pressure stimuli for 3 s to the 
right index fingertip of healthy participants and acquired functional magnetic resonance imaging (fMRI) data using 
a 3T MRI system. 

Results: DCM analysis revealed intra-hemispheric effective connectivity between the contralateral SI (cSI) and Sll 
(cSII) characterized by both parallel (signal inputs to both cSI and cSII) and serial (signal transmission from cSI to 
cSII) pathways during pressure stimulation. DCM analysis also revealed inter-hemispheric effective connectivity 
among cSI, cSII, and the ipsilateral Sll (iSII) characterized by serial (from cSI to cSII) and Sll-level (from cSII to iSII) 
pathways during pressure stimulation. 

Conclusions: Our results support a hierarchical somatosensory network that underlies processing of low-frequency 
tactile information. The network consists of parallel inputs to both cSI and cSII (intra-hemispheric), followed by serial 
pathways from cSI to cSII (intra-hemispheric) and from cSII to iSII (inter-hemispheric). Importantly, our results suggest 
that both serial and parallel processing take place in tactile information processing of static mechanical stimuli as well 
as highlighting the contribution of callosal transfer to bilateral neuronal interactions in Sll. 



Background 

Four types of mechanosensitive afferents mediate the sense 
of touch in the human skin including slow-adapting 
type I (SA-I) afferents for low-frequency (static) stimuli, 
slow- adapting type II (SA-II) afferents for skin stretching, 
fast-adapting type I (FA-I) afferents for flutter, and fast- 
adapting type II (FA-II) afferents for high-frequency 
(vibratory) stimuli [1,2], During mechanical stimulation. 
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they project sensory signals to the somatosensory cortical 
regions for tactile perception [3]. Among them, the som- 
atosensory projection from SA-I afferents is the most ef- 
fective and reliable for object grasping and manipulation 
[4-6] with its characteristic spatial responses [7-9]. 

It is of great interest to understand how tactile informa- 
tion is processed over somatosensory cortical networks. 
Hierarchical organization of tactile information processing 
in the primary (SI) and secondary (SII) somatosensory cor- 
tices has been documented in many anatomical [10-12] 
and neuroimaging studies [13-17]. However, there is on- 
going debate concerning whether tactile information is 
processed in serial (relay of sensory signals from SI to SII) 
or in parallel (relay of sensory signals to both SI and SII). 
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Recently, dynamic causal modeling (DCM) with func- 
tional magnetic resonance imaging (fMRI) data has 
been proposed to address this issue. DCM treats the 
brain as a dynamic system to identify effective connect- 
ivity among brain regions, i.e. which regions cause ac- 
tivity in target regions based on model estimation and 
regional coupling parameters. In DCM, serial process- 
ing models for somatosensory cortical networks 
hypothesize sequential transduction of sensory inputs 
from SI to SII [15], whereas parallel processing models 
hypothesize bifurcated transduction of sensory inputs 
to both SI and SII [16]. A recent fMRI study by Liang 
et al. suggested that parallel processing models may 
better explain effective connectivity in SI and SII for 
electrical and thermal stimuli [16]. In contrast, another 
fMRI study by Kalberlah et al. suggested that serial pro- 
cessing models may better explain effective connectiv- 
ity in SI and SII for vibrotactile stimuli [15]. Thus, 
these reports suggest that a hypothesis for sensory sig- 
nal transduction between SI and SII should be evalu- 
ated depending upon the type of tactile stimuli. In this 
study, we focused on tactile information processing 
models for static (e.g. pressure) stimuli. In addition to 
previous studies that only examined intra-hemispheric 
networks, we further aimed to investigate inter-hemi- 
spheric networks across SI and the bilateral SII. To our 
knowledge, no human fMRI study has investigated ef- 
fective connectivity across hemispheres for any mech- 
anical stimulus. 

In the present study, we addressed two questions re- 
garding somatosensory networks associated with SA-I 
afferents: first, we addressed how intra-hemispheric ef- 
fective connectivity is formed in the contralateral SI 
(cSI) and SII (cSII) for tactile information processing of 
pressure stimuli. We employed DCM to clarify whether 
pressure stimuli are processed in serial or in parallel 
with similar hypotheses to previous reports [15,16]: 
(1) the serial processing model hypothesis highlighting 
sequential inputs from cSI to cSII and (2) the parallel 
processing model hypothesis highlighting two-way in- 
puts to cSI and cSII. Second, we addressed how inter- 
hemispheric effective connectivity is formed across cSI, 
cSII, and the ipsilateral SII (iSII) for tactile information 
processing of pressure stimuli. Consequently, we per- 
formed a second DCM analysis to assess three possible 
models: (1) the first model hypothesized causal activity 
(an information flow) from cSI to iSII; (2) the second 
model hypothesized causal activity from cSII to iSII; 
and (3) the third model hypothesized causal activity 
from cSI to iSII as well as from cSII to iSII. Using the 
human fMRI data recorded from our pressure stimula- 
tion experiment, we evaluated each hypothesis to find 
which sensory signal transduction model most likely 
explained the neural data. 



Methods 

Participants 

Twenty-one healthy volunteers (age, 24.19 ± 2.71 years; 
all right-handed) with no history of neurological disor- 
ders or deficits in tactile processing gave written in- 
formed consent and participated in this study, which 
was approved by the Korea University Institutional Re- 
view Board (KU-IRB-11-46-A-1). 

Pressure stimulation 

A band-type MR-compatible stimulation device built by 
our group (Figure 1) [18] was used to apply a pressure 
stimulus of 5.98 g/mm^ to the right index fingertip. A 
neonatal cuff (Ml 866 A, Philips Healthcare, Best, The 
Netherlands) wrapped around the fingertip was directly 
connected to a rolling pump in a blood pressure monitor 
(BP3AG1, Microlife AG, Widnau, Switzerland) through 
an elastic air-tube with a length of 5 m and a diameter 
of 4 mm. The cuff was controlled by a pressure sensor 
for achieving uniform pressure; the sensor in turn was 
controlled by E-Prime 2.0 software (Psychology Software 
Tools, Inc., Sharpsburg, PA, USA) for configuring the 
length of stimulation. The cuff expanded at the turn-on 
of the pump and pressed the whole ventral surface of 
the fingertip (see Kim et al. for further information con- 
cerning the stimulation device [18]). 

Participants were instructed to lie comfortably on the 
MRI table with eyes closed during the scanning in each 
session, to put earplugs in their ears, and to hold an 
emergency squeeze-ball in their left hands during the en- 
tire scanning session. Before scanning, the cuff was at- 
tached to wrap participants' right index fingertips with 
minimal pressure on the skin. Each participant performed 
four block-designed study sessions. Each session consisted 
of four trials. To avoid potential adaptation due to repeti- 
tive stimulus application, we designed four separate ses- 
sions instead of one session, totaling 16 trials. A single trial 
comprised a 21 s resting period followed by a 3 s stimula- 
tion period. During each stimulation period, a single static 
indentation was applied continuously to the participant's 
right index fingertip. 

Anatomical and functional data acquisition 

Anatomical and functional images were acquired using 
a 3T MRI system (Magnetom TrioTim, Siemens Med- 
ical Systems, Erlangen, Germany) with a standard 32- 
channel head coil. Ti-weighted anatomical images were 
acquired using a 3D magnetization-prepared gradient 
echo (MPRAGE) sequence, with the imaging parameters 
of repetition time (TR) = 1,900 ms, echo time (TE) = 
2.48 ms, flip angle = 9°, field of view (FOV) = 200 mm, 
and voxel size = 0.8 x 0.8 x 1 mm^. T2-weighted functional 
images were acquired using a gradient echo-planar im- 
aging (EPI) sequence, with the imaging parameters of 
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Figure 1 An MR compatible stimulation device. A neonatal cuff (A) wrapped around the right index fingertip (B: front view; C: side view) 

exerted a pressure stimulus of 5.98 g/mm^ by expanding at the turn-on of the pump in a blood pressure monitor through an elastic air-tube. It 

pressed the whole ventral surface of the fingertip for 3 s in each stimulus period. 

J 



TR = 3,000 ms, TE = 30 ms, flip angle = 90°, FOV = 
240 mm, slice thickness = 3 mm, and voxel size = 3 x 
3x3 mm^. 

Statistical analysis 

Functional images were preprocessed using SPM8 
(Wellcome Department of Imaging Neuroscience, UCL, 
London, UK), through a series of steps of slice-timing 
correction, realignment with the rigid-body transform- 
ation matrices, normalization to the Montreal Neuro- 
logical Institute (MNI) template, and smoothing with an 
8 mm full-width-half-maximum (FWHM) isotropic 
Gaussian kernel The mean EPI image of each individual 
subject was directly warped into the standard EPI tem- 
plate in SPM8 during the normalization step. Then, the 
conventional general linear model (GLM) in SPM8 per- 
formed statistical analyses on blood oxygenation level- 
dependent (BOLD) signals with a canonical hemodynamic 
response function and its time and dispersion derivatives. 
A 128 s high-pass filter removed physiological artifacts 
in the BOLD signals. A full factorial design based on a 
random effects model performed a group analysis for 
the inference of statistically significant cortical activa- 
tion. Cluster-level F-statistics {p < 0.05 with a family- 
wise-error (FWE) correction, with a minimum threshold 
(k) of 5 voxels for significant clusters) produced group- 
level statistical parametric maps (SPMs) representing 
significant voxel clusters. The automated anatomical la- 
beling toolbox [19] determined anatomical cluster labels 
of the activation regions in the SPMs. 

DCM analysis 

We used the Anatomy toolbox [20] to generate anatom- 
ical masks for three seed regions of interests (ROIs), in- 
cluding cSI (Brodmann area (BA) 3a, 3b, 1, and 2) [21,22], 
cSII, and iSII (parietal operculum (OP) 1, 2, 3, and 4) 
[23,24] for DCM analysis (Figure 2A). DCM infers effect- 
ive connectivity by estimating parameters of regional 
coupling using a Bayesian framework in dynamic systems 



of neuronal populations, which are unobservable directly 
from BOLD signals. DCM is a hypothesis -driven approach, 
and finds optimal model parameters at the neuronal level 
to make BOLD signals generated with predefined hy- 
potheses as close to observed BOLD signals as possible. 
We used a bilinear state equation with three components: 
(1) experimental (driving) inputs perturbing brain states, 
i.e., in our case, sensory signals directly projecting to the 
cortex; (2) intrinsic connectivity in the absence of expe- 
rimental perturbations; and (3) changes (modulations) of 
the intrinsic connectivity induced by experimentally ma- 
nipulated inputs, i.e. changes in regional couplings by sen- 
sory inputs, which provided information concerning how 
much activation in source regions receiving direct inputs 
caused an increase/decrease in activation in target regions 
per unit of time [25] . The resulting modulations were the 
components of interest in this study and were used to 
model the flow of tactile information among the somato- 
sensory cortical regions. Bayesian model selection (BMS) 
in DCM was achieved by a free energy approximation to 
the log evidence of each model in terms of model fit and 
complexity. BMS determined which model was the most 
preferred among defined models [26,27]. 

To model intra- and inter-hemispheric effective connect- 
ivity in this study, we used DCM instead of other methods 
for effective connectivity such as Granger causality [28] 
and structural equation modeling [29] because DCM has 
some advantages that (1) it provides direct modeling of ef- 
fective connectivity between brain regions at the hidden 
neural level [25]; (2) it explains the relationship between 
neural activity and BOLD hemodynamic responses [25,26]; 
(3) it is less liable to be affected by the variability of 
hemodynamic response functions [30,31]; and (4) it pro- 
vides modeling of effective connectivity caused by experi- 
mental perturbations [25,31]. 

Three seed ROIs were selected in each participant from 
individual analyses using the anatomical masks with 
cluster-level F-statistics (uncorrected 0.001, k>5). 
Among 21 participants, six participants who showed 
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Figure 2 ROIs and model families for DCM analysis. (A) Three seed ROIs are shown in cSI, cSII, and iSII (group analysis with 21 participants, 
F-test, P(FWE) < 0.05, k>5, bar: F-statistics) (B) Two intra-hemispheric model families were constructed to find out whether pressure stimuli were 
processed in serial (left: family A) or parallel (right: family B) in cSI and cSII. With full intrinsic connectivity (thin arrows), we defined at least one 
and at most four modulations (empty circles) in the intrinsic connections. A key difference between family A and family B was that one driving 
input (a thick arrow) was applied to only cSI in family A, whereas two driving inputs (two thick arrows) were applied to both cSI and cSII. (C) Three 
inter-hemispheric model families were constructed to find out how pressure stimuli were delivered within cSI, cSII, and iSII: from cSI to iSII (left), from 
cSII to iSII (middle), or both (right). With full intrinsic connectivity, we defined a fixed modulation (filled circles) from cSI to iSII but no modulation from 
cSII to iSII (left), a fixed modulation from cSII to iSII but no modulation from cSI to iSII (middle), and two fixed modulation from cSI to iSII and from cSII 
to iSII (right). At least one modulation (empty circles) was applied to the intrinsic connections. Two driving inputs (two thick arrows) were applied to 
both cSI and cSII. 



no significant cluster for at least one of the three seed 
ROIs were excluded. Then, in each of the remaining 15 
participants, we defined volumes of interests (VOIs) as 
spheres of 6 mm radius centered on the most signifi- 
cant peaks of individual ROIs. From the BOLD signals 
of all voxels in each VOI, the first eigenvariate was ex- 
tracted as a representative time-series. Table 1 shows 
the average MNI coordinates and the average numbers 
of voxels across participants in cSI, cSII, and iSII. 

We performed two DCM analyses. First, we con- 
structed two intra-hemispheric model families corre- 
sponding to serial processing (family A) or parallel 
processing (family B) between cSI and cSII. We defined 
full intrinsic connectivity due to the fact that both cSI 
and cSII were fully anatomically connected [32,33]. It 
led to four intrinsic connections: self-connections in 
each of cSI and cSII, a connection from cSI to cSII, and 
a connection from cSII to cSI. We presumed that at least 
one and at most four connections could be simultan- 
eously modulated, as tactile sensory inputs would affect 
the dynamics of at least one intrinsic connection. This 



resulted in 15 (2^ - 1) models of modulations excluding 
a case with no modulation. Thus, there were a total of 
30 models including 15 in family A and 15 in family B. 
A key difference between model families was that one 
driving input (a pressure stimulus) was applied only to 
cSI in family A, whereas two driving inputs were applied 
both to cSI and cSII in family B. In the present study, 
we examined a total of 450 models (combinations of 2 
famiUes, 15 models, and 15 participants) (Figure 2B). 



Table 1 The average MNI coordinates and the average 
number of voxels in VOIs (15 participants) 



vols 


Average MNI coordinates (mm) 


Average number of voxels 




X 


y 


z 




cSI 


-51 (±1) 


-20 (±2) 


49 (±2) 


25 (±2) 


cSII 


-50 (±2) 


-15 (±1) 


17 (±1) 


27 (±1) 


iSII 


52 (±2) 


-13 (±1) 


16 (±1) 


23 (±2) 



cSI: contralateral primary somatosensory cortex; cSII: contralateral secondary 
somatosensory cortex; iSII: ipsilateral secondary somatosensory cortex; 
numbers in parentheses: standard error. 
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Next, we constructed three inter-hemispheric model 
famiUes to infer how the information of pressure stimuU 
was conveyed among cSI, cSII, and iSIL We again defined 
full intrinsic connectivity between all pairs of three VOIs 
due to the fact that cSI, cSII, and iSII were fully anato- 
mically connected [34,35]. This resulted in 9 connec- 
tions (3 self-connections and 6 bidirectional connections 
between three VOIs) with 512 (= 2^) modulation models. 
We applied two driving inputs to both cSI and cSII. To 
explain tactile information flows over inter-hemispheric 
connections, we built three model families. Each family 
consisted of a set of models with specific modulations 
being fixed: those fixed modulations represented pre- 
determined information pathways in accordance with 
our three hypotheses. We hypothesized that iSII re- 
ceived tactile information from cSI in family C (cSI^ 
iSII), from cSII in family D (cSII iSII), and from both 
cSI and cSII in family E (both cSI^iSII and cSII ^ 
iSII). Composition of modulation models in each family 
was determined as follows. 

In family C, we first determined a fixed modulation in 
cSI^iSII and no modulation in cSII^iSII, which left 
128 (= 2^) out of 512 models to be configured. Similarly, 
there were also 128 configurable models by determining 
a fixed modulation in cSII iSII and no modulation in 
cSI iSII in family D, or by determining fixed modula- 
tions in both cSI iSII and cSII iSII in family E. 
Then, among 128 models in each family, we excluded 
those models that had modulation in neither cSI^cSII 
nor cSI <— cSII as we avoided a case when there was no 
intra-hemispheric connectivity but inter-hemispheric con- 
nectivity. In other words, we considered integrating inter- 
hemispheric connectivity models with intra-hemispheric 
models by assuming that inter-hemispheric connectivity 
should involve intra-hemispheric connectivity. This exclu- 
sion procedure removed 32 models (= 2^) from 128, 
resulting in 96 models for each family. Thus, in our inter- 
hemispheric DCM analysis, we examined a total of 4320 
models (combinations of 3 families, 96 models, and 15 
participants) per session (Figure 2C). 

To assess the fitness of each model to our experimental 
data, we performed BMS based on a random effect infe- 
rence assuming that model structures could vary across 
participants [26,27]. BMS determined the best model and 
the best model family by computing exceedance proba- 
bilities of all models and model families. The exceedance 
probability inferred the probability that a specific 
model (or family) described the data better than any 
other model (or family) being compared [36]. Finally, 
the influence of modulation in the best models were 
evaluated using ^-tests to ascertain whether modula- 
tory parameters were statistically significant across par- 
ticipants with the null hypothesis being that differences 
were equal to zero. 



Results 

As we used simple and distinct pressure stimuli with 
forces sufficiently higher than the absolute sensitivity 
threshold for the fingertips [37], all participants clearly 
felt the induced static indentation on their ventral sur- 
faces of the right index fingertips for 3 s and confirmed 
this sensation in post-hoc interviews. Table 2 shows de- 
tailed information concerning significant clusters of acti- 
vation during pressure stimulation (F-test, P(ywe) < 0-05, 
k>5). These results highlighted that cortical regions 
known to be related to tactile perception were activated 
in response to pressure stimulation during our experi- 
ment (Figure 3). 

The DCM analysis for intra-hemispheric effective con- 
nectivity revealed that family B (parallel processing) was 
preferred to family A (serial processing). BMS resulted in 
exceedance probabilities for family A and B of 13.28% and 
86.72%, respectively (averaged across sessions). Among 15 
single models in family B, the best model showed an ex- 
ceedance probability of 30.27% (averaged across sessions). 
This model contained one modulation in a forward connec- 
tion from cSI to cSII (a serial pathway: cSI^cSII). We 
confirmed the statistical significance of the modulatory pa- 
rameters from cSI to cSII (^-test, p < 0.0005) across parti- 
cipants and sessions using the best single model (Figure 4). 

DCM analysis for inter- hemispheric effective connec- 
tivity revealed that family D (cSII^iSII) was the most 
preferred model family. BMS resulted in exceedance pro- 
babilities for family C, D, and E of 31.17%, 58.30%, and 
10.53%, respectively (averaged across sessions). Among 96 
single models in family D, we found the best model with 
an exceedance probability of 7.77% (averaged across ses- 
sions). This model contained two modulations, one in a 
forward connection from cSI to cSII (a serial pathway: 
cSI cSII) and the other in a forward connection from 
cSII to iSII (a serial pathway: cSII iSII). We tested the 
statistical significances of the modulatory parameters in 
the best single model from cSI to cSII (^-test, p < 0.0005) 
and those from cSII to iSII (^-test, p < 0.05) across par- 
ticipants and sessions (Figure 5). 

Discussion 

Intra-hemispheric effective connectivity 

Considering the ongoing debate concerning serial and 
parallel processing in tactile perception, the results ob- 
tained in this study provide support for both hypotheses. 

Prior evidence for the parallel processing hypothesis 
includes the anatomical thalamic neuronal projection 
from the ventroposterior thalamus (VPL) to both SI and 
SII in cats [38] and monkeys [39], robust activation of 
SII unaffected by inactivation of SI during tactile stimu- 
lation in cats [40], and simultaneous response onsets of 
activities in SI and SII to laser-induced nociceptive stim- 
uli in a human magnetoencephalography (MEG) study 
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Table 2 Significantly activated clusters during pressure stimulation (group analysis with 21 participants, F-test, 
P(FWE) < 0.05, /c> 5) 



Anatomical labels 




X 


MNI coordinates (mm) 

y 


z 


Voxels 


F 


Z 


P(FWE) 


Insula 


R 


42 


3 


9 


116 


38.41 


7.43 


0.000 


Insula 


R 


39 


-3 


-3 




23.81 


6.20 


0.000 


Postcentral gyrus 


L 


-51 


-21 


18 


527 


32.57 


7.00 


0.000 


Insula 


L 


-39 


-3 


9 




30.26 


6.81 


0.000 


Inferior parietal lobule 


L 


-57 


-21 


45 




23.71 


6.19 


0.000 


Rolandic operculum 


R 


54 


-21 


21 


220 


26.79 


6.50 


0.000 


Postcentral gyrus 


R 


60 


-18 


33 




19.38 


5.68 


0.000 


Postcentral gyrus 


R 


51 


-21 


45 




14.60 


4.99 


0.017 


Precentral gyrus 


R 


30 


-24 


60 


207 


23.12 


6.13 


0.000 


Precentral gyrus 


R 


24 


-15 


72 




20.81 


5.86 


0.000 


Paracentral lobule 


R 


12 


-27 


72 




18.51 


5.57 


0.001 


Median cingulate 


L 


-6 


9 


36 


6 


16.85 


5.34 


0.003 


Insula 


R 


33 


27 


3 


22 


16.61 


5.30 


0.004 


Insula 


L 


-30 


21 


6 


8 


15.61 


5.15 


0.008 


Paracentral lobule 


L 


-9 


-27 


66 


11 


14.52 


4.98 


0.018 



[41]. Conversely, prior evidence for the serial processing 
hypothesis includes the anatomical serial connections of 
neurons from SI to SII in monkeys [11,12], sequential 
response onsets from SI to SII with temporal differences 
between the regional onsets of < 100 ms to electrical 
stimuli in a human MEG study [13,42,43], and three se- 
quential information channels from the thalamus to cSI, 
from cSI to cSII, and from cSII to iSII based on response 
onsets in cSI, cSII, and iSII during electrical stimulation 
in a human electroencephalography (EEG) study [14]. 

A review of somatosensory evoked potentials (SEPs) 
highlighted that SII directly received nociceptive inputs 
from the thalamus, whereas it mainly received tactile 
inputs from SI [44]. This report hinted at the possibility 
of the coexistence of the two hypotheses. Additionally, 
two human fMRI studies using DCM supported each 



hypothesis; one study reported that nociceptive (laser- 
induced heat) and non-nociceptive (electrical pulses) 
stimuli were directly delivered to both cSI and cSII 
through parallel processing [16], whereas the other 
study reported that vibrotactile stimuli were transmitted 
from cSI to cSII through serial processing [15]. Based on 
these previous studies, we anticipated the possibility of the 
coexistence of serial and parallel characteristics in tactile 
information processing for non-nociceptive mechanical 
stimuli. Our first DCM results are consistent with a study 
of Liang et al. [16] in terms of parallel inputs of sensory sig- 
nals to cSI and cSII (family B), and also in agreement with 
a study of Kalberlah et al [15] in terms of the serial trans- 
mission of sensory signals from cSI to cSII (the best model 
in family B contained a single modulation connection from 
cSI to cSII). However, our results do not perfectly connect 




Figure 3 Cortical activation during pressure stimulation. GLM analysis revealed significantly activated clusters in the insula, postcentral gyrus, 
rolandic operculum, supramarginal gyrus, precentral gyrus, median cingulate, and paracentral lobule (group analysis with 21 participants, F-test, 
P(FWE) < 0.05, k> 5, bar: F-statistics, A: sagittal view, B: coronal view, C: transverse view). The left side at B and the upper side at C indicate the 
contralateral hemisphere (contra.); the right side at B and the bottom side at C indicate the ipsilateral hemisphere (ipsi.); and activation at A 
indicates all suprathreshold voxels from both the contralateral and ipsilateral hemispheres. 

V J 



Chung et al. BMC Neuroscience 2014, 15:43 
httpy/www.biomedcentral.com/l 471-2202/1 5/43 



Page 7 of 10 




Family A Family B 

30 models in total 



Figure 4 DCM results for intra-hemispheric effective connectivity. 

Family B (parallel processing, exceedance probability of 86.72%) was 

preferred to family A (serial processing, exceedance probability of 

13.28%) (top left). The best single model (model 17) with an 

exceedance probability of 30.27% among 15 models in family B was 

selected (bottom). The best single model contained a modulation in 

the connection from cSI to cSII, implying a serial pathway. Numbers 

along the intrinsic connections indicate the mean parameters of 

intrinsic connectivity, and numbers in parentheses indicate the mean 

parameter of modulation from cSI to cSII across participants and 

sessions of the best single model (top right). The increase in the 

activity of cSII corresponds to 43% of the activity of cSI by pressure 

stimulation (0.43 in the parentheses from cSI to cSII). 
\ J 



the two previous DCM studies as Liang et al considered 
models with driving inputs of electrical stimuli to the thal- 
amus and Kalberlah et al used vibrotactile stimuli. Thus, 
our study may reveal intra-hemispheric effective connectiv- 
ity for another type of sensory signals from SA-I afferents 
characterized by: (1) the parallel processing of a bifurcated 
sensory input to cSI and cSII; and (2) the serial processing 
of sequential signal transduction from cSI to cSII. In terms 
of the afferent-dependency of intra-hemispheric effective 
connectivity based on previous and our DCM analyses, we 
conjecture that very high-frequency electrical stimulation 
(delivered by FA-II) results in parallel processing (from 
Liang et al. 2011 [16]), high-frequency vibrotactile stimula- 
tion (delivered by FA-I) results in serial processing (from 
Kalberlah et al. 2013 [15]), whereas very low-frequency 
vibrotactile stimulation (delivered by SA-I) recruits both 
processing types (from our study). 




Family C Family D Family E 

288 models in total 



Figure 5 DCM results for inter-hemispheric effective connectivity. 

Family D (Sll-level information flow, exceedance probability of 

58.30%) was the most preferred model family, compared with 

family C (exceedance probability of 31.17%) and family E 

(exceedance probability of 10.53%) (top left). The best single model 

(model 2) with an exceedance probability of 7.77% among 96 

models in family D was selected (bottom). The best single model 

contained two modulations in the connections from cSI to cSII (the 

serial pathway shown in intra-hemispheric effective connectivity) and 

from cSII to iSII. Numbers along the intrinsic connections indicate the 

mean parameters of intrinsic connectivity and numbers in parentheses 

indicate the mean parameters of modulations from cSI to cSII and from 

cSII to iSII across participants and sessions of the best single model 

(top right). The increase in the activity of cSII corresponds to 32% of 

the activity of cSI (0.32 in the parentheses from cSI to cSII) and the 

increase in the activity of iSII corresponds to 16% of the activity of cSII 

(0.16 in the parentheses from cSII to iSII) by pressure stimulation. 
\ J 

Inter-hemispheric effective connectivity 

Anatomical studies in cats [45] and monkeys [46] have 
shown the presence of callosal projection neurons in SI 
and SII and postulated inter-hemispheric transfer of 
somatosensory information for body representation. An 
animal study in monkeys revealed denser callosal con- 
nections from BA 3b to 2 implying an inter-hemispheric 
pathway originating from cSI [47]. A human MEG study 
suggested almost simultaneous delivery of sensory sig- 
nals of electric stimuli to cSII and iSII in a latency < 
4 ms [48]. However, how tactile information is delivered 
from one side to the other remains controversial. It has 
been reported that SII has sensory neurons with bilateral 
receptive fields and dense callosal fibers [49]. Thus, a num- 
ber of studies have suggested that the corpus callosum 
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supports bilateral receptive fields in SII neurons for inter- 
hemispheric tactile information transfer. Animal studies 
have demonstrated decreases in the proportion of SII neu- 
rons with bilateral receptive fields in callosotomized cats 
with normal cats [50,51]. Human fMRI studies have dem- 
onstrated involvement of SII and the posterior parietal cor- 
tex (PPC) in inter-hemispheric tactile information transfer, 
by showing absence of activation responding to tactile 
stimuli in the ipsilateral SII and PPC in fully or partially 
callosotomized patients [52]. 

Based on these studies, we derived two hypotheses of 
inter-hemispheric pathways from cSI to iSII (originating 
from cSI) or from cSII to iSII (within SII). A recent SEP 
study demonstrated time-varying source connectivity dur- 
ing electrical stimulation such that tactile information flo- 
wed from cSI to iSII followed by information flow from 
cSII to iSII with a short latency of approximately 15 ms 
[14]. Consequently, we advanced another hypothesis of an 
inter-hemispheric pathway combining the other two hy- 
potheses, both from cSI to iSII and from cSII to iSII. 
Among the three hypotheses, our DCM analysis results 
support the second hypothesis representing an Sll-level 
inter-hemispheric pathway from cSII to iSII (family D). 
Additionally, the best model in family D supported serial 
processing as shown in the first DCM results. 

The Sll-level pathway for tactile information delivery we 
found may be attributed to the roles of SII for high-level 
tactile perception. Animal studies in monkeys reported the 
roles of SII for tactile discrimination and learning by in- 
vestigating impaired task performance of monkeys after 
removal of the bilateral SII [53-55]. A recent human 
MEG-fMRI study implicated a callosal interconnection 
of the bilateral SII in bimanual tactile exploration of ob- 
jects, revealing a significant relationship between task 
performance and inter-hemispheric inhibition within 
SII for encoding or comparing spatial features of objects 
[49]. Another human fMRI study suggested more com- 
plex bilateral receptive fields in SII than in SI for per- 
ceiving higher order features of tactile stimuli based on 
the investigation of distinct bilateral activation in SII to 
mechanical stimuli on fingers [56]. In addition, the role 
of the bilateral interaction within SII was explained by 
sensory-motor integration for precise control of move- 
ment executions because of multiple reciprocal connec- 
tions between SII and other cortical regions (e.g. motor 
areas) [51]. 

Based on our results of two DCM analyses, we suggest 
that tactile information delivered through SA-I afferents 
by pressure stimulation is processed not only in parallel 
with a bifurcated input to both cSI and cSII, but also in 
serial with sequential sensory transmission from cSI to 
cSII intra-hemispherically and from cSII to iSII inter- 
hemispherically. These sensory information pathways 
are consistent with previous results that include sensory 



projections from thalamus to both cSI and cSII, hier- 
archical sensory information processing from cSI to cSII, 
and bilateral interactions by bilateral receptive fields 
within SII. Our results of inter-hemispheric connectivity 
underline this bilateral interaction of SII for tactile infor- 
mation processing by showing that pressure stimulation 
information from cSI more likely flows through cSII-iSII 
connections, not directly to iSII. Therefore, our intra- 
and inter-hemispheric effective connectivity model can 
be considered as an integrated model, which combined 
the three architectural factors explained above, can be 
used to understand low- and high-level tactile informa- 
tion processing of pressure stimulation. 

Limitations and future work 

In our investigation of hierarchical somatosensory net- 
works, we ruled out any temporal variation because the 
state equation in DCM did not consider inter-regional 
conduction delays [26]. We explained here only inter- 
regional causal activities implying information flow, not 
any timing issue, e.g. activation onset in cSI, cSI, and 
iSII. Thus, further studies are needed to investigate tem- 
poral characteristics to corroborate our effective con- 
nectivity models in view of temporal information flow 
during pressure stimulation. 

In addition, our DCM analysis needs to be extended 
to other mechanical stimuU (flutter or high-frequency 
vibration) from FA afferents, electrical stimuU, and even 
nociceptive stimuli to evaluate the modality-specific 
consistency of our intra- and inter-hemispheric effective 
connectivity models. Until now, it has been reported that 
mechanical stimuli from FA afferents were processed in 
serial [15] and nociceptive stimuli were processed in paral- 
lel [16,41,57,58]. However, both serial [13,14,59] and paral- 
lel [16,39] processing modes were reported in separate 
studies for electrical stimulation. Our study suggests that 
mechanical stimuli, limited to signals from SA-I afferents 
(pressure stimuli), were processed both in serial and paral- 
lel. Hence, further evaluation of our models is required 
with other sensory stimuli from different afferents. 

In terms of higher-level tactile processing, a recent 
neuroimaging study reported the existence of dual som- 
atosensory pathways for the perception of texture (from 
SI to SII) and location (directly to SII) [60]. It therefore 
would be interesting to use such high-level tactile stimuli 
(e.g. texture or location) in addition to low-level tactile 
stimuli (static indentation as done here) to validate the 
effectiveness of our models for information processing 
of low-level tactile perception in the context of higher- 
level tasks. 

Finally, causal relationships between cortical regions 
estimated from the fMRI data should also be examined 
using other assessment methods, in particular model- 
free methods for connectivity such as the information- 
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theoretic transfer entropy method [61,62] to justify 
model-free measures of information flow. 

Conclusions 

In the present study, we investigated somatosensory net- 
works based on effective connectivity for information pro- 
cessing of tactile signals from SA-I afferents. Our first 
DCM analysis revealed that tactile signals were processed 
through intra-hemispheric effective connectivity charac- 
terized by parallel (sensory inputs to cSI and cSII) and 
serial (signal transmission from cSI to cSII) pathways, sup- 
porting both serial and parallel processing of tactile infor- 
mation. Our second DCM analysis revealed that tactile 
signals were transmitted serially from cSI, through cSII, to 
iSII over inter-hemispheric connections. Consequently, we 
postulate a tactile information pathway of pressure stimuli 
in cSI, cSII, and iSII with three components: (1) parallel 
processing with a bifurcated input to both cSI and cSII; 
(2) a serial pathway from cSI to cSII; and (3) an Sll-level 
pathway from cSII to iSII. 
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